First-principles LDA+L7 and GGA+C/ study of plutonium oxides 
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The electronic structure and properties of Pu02 and PU2O3 have been studied from first principles 
by the all-electron projector-augmented-wave (PAW) method. The local density approximation 
(LDA)+(7 and the generalized gradient approximation (GGA)+(/ formalism have been used to 
account for the strong on-site Coulomb repulsion among the localized Pu 5/ electrons. We discuss 
how the properties of PUO2 and PU2O3 are affected by the choice of U as well as the choice of 
exchange-correlation potential. Also, oxidation reaction of PU2O3, leading to formation of Pu02, and 
its dependence on U and exchange-correlation potential have been studied. Our results show that by 
choosing an appropriate U it is promising to correctly and consistently describe structural, electronic, 
and thermodynamic properties of PUO2 and PU2O3, which enables it possible the modeling of redox 
process involving Pu-based materials. 



I. INTRODUCTION 

Plutonium dioxide (PUO2) and sesquioxide (PU2O3) 
are the only observed stoichiometric compounds formed 
at the surface of the metallic plutonium when exposed 
to dry air— (nonstoichiometric oxide may form by reac- 
tion of dioxide with water 2 ). From this sense, pluto- 
nium corrosion and oxidation are often treated as equiv- 
alent topic. The plutonium corrosion plays a key role in 
considering the nuclear stockpile and storage of surplus 
plutonium. Therefore, a thorough understanding of the 
physical and chemical properties of the plutonium oxide 
is always needed. 

From basic point of view, it can be visualized that 
many physical and chemical properties of the plutonium 
oxide are closely related to the quantum process of lo- 
calization and delocalization for Pu 5/ electrons. Model- 
ing of the electron localization, and thus any redox pro- 
cess involving plutonium, is a complex task. Conven- 
tional density functional schemes that apply the local 
density approximation (LDA) or the generalized gradi- 
ent approximation (GGA) underestimate the strong on- 
site Coulomb repulsion of the Pu 5/ electrons and conse- 
quently fail to capture the correlation-driven localization. 
Therefore, the 5/ electrons in elemental Pu, as well as in 
Pu compounds, require special attention. One promis- 
ing way to improve contemporary LDA and GGA ap- 
proaches is to modify the intra-atomic Coulomb interac- 
tion through the so-called LDA+[/ or GGA+C/ approach, 
in which the underestimation of the intraband Coulomb 
interaction is corrected by the Hubbard U parameter^! 5 -. 
This method has been used to discuss the equilibrium 
lattice parameter of bulk Pu in Ref.-i£. The choice of 
U is, however, not unambiguous and it is not trivial to 
determine its value a priori, though there are attempts 
to extract it from standard first-principles calculations. 
Hence, U is often fitted to reproduce a certain set of ex- 
perimental data such as band gaps and structural prop- 
erties. 

In this paper we use the LDA+J7 and GGA+J7 schemes 
due to Dudarev et alM- to calculate the lattice parame- 
ters, electronic structure, and thermodynamic properties 



of PUO2 and PU2O3. We discuss how these properties 
are affected by the choice of U as well as the choice 
of exchange-correlation potential, i.e., the LDA or the 
GGA, and how redox processes occurred in plutonium ox- 
ide can be explored in the LDA+?7 and GGA+C/ formal- 
ism. In addition, we notice that recently there have oc- 
curred a few experimenta l 10 ! 11 ' 12 and theoretica l 11 ' 13 ' 14 ' 15 
studies of the electronic structures of plutonium ox- 
ides. In this paper we have compared our calculated 
LDA/GGA+C/ results with those reports. Our results 
show that while the pure LDA/GGA (without U correc- 
tion) calculations fail to describe the ground-state be- 
haviors of the plutonium oxides, such as the insulat- 
ing nature, the magnetic configuration, and the 5/ band 
gap, the present LDA/GGA+C7 approaches with tunable 
Coulomb parameters can effectively remedy those fail- 
ures and the consequent results fit well in the attainable 
experimental dat a 10 ' 11 ' 12 . 

This paper is organized as follows. The details of our 
calculations are described in Sec. II and in Sec. Ill we 
present and discuss the results. In Sec. IV, we summarize 
our findings. 



II. METHODOLOGY OF THE CALCULATION 



The calculations were performed using the projector- 
augmented wave (PAW) method of Blochl 16 , as im- 
plemented in the ab initio total-energy and molecular- 
dynamics program VASP (Vienna ab initio simulation 
program) 17 ' 18 ' 19 ' 20 . PAW is an all-electron method that 
combines the accuracy of augmented-plane-wave meth- 
ods with the efficiency of the pseudopotential approach. 
The PAW method is implemented in VASP with the 
frozen-core approximation. For the plane- wave set, a cut- 
off energy of 400 eV was used. The plutonium 6s, 6p, 7s, 
and 5/, and the oxygen 2s and 2p electrons were treated 
as valence electrons. The strong on-site Coulomb repul- 
sion amongst the localized Pu 5/ electrons are accounted 
for by using the formalism formulated by Dudarev et aZ3. 
In this scheme the total LDA (GGA) energy functional 
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FIG. 1: (a) Unit cell of PUO2 containing 3 atoms. The black 
spheres are Pu atoms, the gray ones are oxygens, (b) Unit 
cell of PU2O3 containing five atoms. 

is of the form 

^LDA(GGA) + (7 = -^LDA(GGA) H ^ ^2 ^P" _ Tr ' P" 

. . (!) 

where p G is the density matrix of / states, and U and 
J are the spherically averaged screened Coulomb energy 
and the exchange energy, respectively. In this paper the 
Coulomb U is treated as a variable, while the exchange 
energy is set to be a constant J=0.75 eV. This value 
of J is in the ball park of the commonly accepted one 
for Pu & 21 ' 22 i 23 . Since only the difference between U and 
J is significant 9 , thus we will henceforth label them as 
one single parameter, for simplicity labeled as U, while 
keeping in mind that the non-zero J has been used during 
calculations. 

The exchange and correlation effects were treated 
in both the local density approximation and the gen- 
eralized gradient approximation 2 ^. We studied PuC>2 
in its ground-state fluorite structure (Fm3m) and the 
sesquioxide PU2O3 in the hexagonal /3-type structure 
(P3ml). For Pu02 we used a 11 x 11 x 11 Monkhorst- 
Pack fc-point mesh 2 ^ (56 irreducible k points) and for 
PU2O3 we used a 9 x 9 x 6 grid (57 irreducible k points). 
The electronic density of states (DOS) was obtained with 
15 x 15 x 15 (120 irreducible k points) and 11x11x9 grid 
(120 irreducible k points) fc-point meshes, respectively. 
The Brillouin-zone integration was performed using the 
modified tetrahedron method of Blochl 26 . In order to 
study the reaction energy it is necessary to calculate the 
energy of an oxygen molecule (Eo 2 )- The effect of spin 
polarization was included in calculating Eo 2 ■ 

III. RESULTS AND DISCUSSION 

A. Atomic and electronic structure of PuC>2 

Plutonium dioxide crystallizes in a CaF2-like ionic 
structure [Fig. 1(a)] with the plutonium and oxygen 
atoms forming face-centered and simple cubic sublattices, 



respectively. In this arrangement each plutonium atom 
is located at the center of an oxygen cube, and for ev- 
ery four such cubes there is an empty one. In the ionic 
limit, formal charge for plutonium in PuC>2 is +4, cor- 
responding to formal population of f 4 . This leads to 
local 5=2 plutonium moment, which can couple with 
other sites in either a ferromagnetic (FM) or antifer- 
romagnetic (AFM) manner. PuC>2 is known to be an 
insulator— and some scattered experimental data 2 ^ sup- 
port the ground state of PuC-2 to be an AFM phase. 
In the present LDA/GGA+t/ approaches, we have con- 
sidered the FM, AFM, and nonmagnetic phases for each 
choice of the value of U and then determined the ground- 
state phase by a subsequent total-energy comparison of 
these three phases. For PUO2, at U=0, the ground state 
is a FM metal, which is in contrast to experiment. By 
increasing the amplitude of U, our LDA/GGA+[7 cal- 
culations correctly predicted an AFM insulating ground 
state. The turning value of U for this FM-AFM energy 
transition of the ground state is of ~1.5 eV. In the dis- 
' cussion that follows, we therefore confined our report to 
the AFM solution for the PuC>2. A thorough description 
of the magnetic structure of plutonium oxides is beyond 
our intention in this paper, and we would like to leave it 
for the future studies. 

The experimentally determined lattice parameter of 
Pu0 2 is a =5.396 Aat 25°CS. Here the calculated a 
and bulk modulus Bq of Pu02 were obtained from the 
corresponding energy minimization at constant volumes 
and by fitting a Murnaghan equation of stated to the re- 
sulting energy- volume data, respectively. The results as a 
function of U within the LDA and the GGA schemes are 
shown in Fig. 2(a) and (b) for ao and Bq, respectively. 
For comparison, the experimental values of and Btp- 
are also shown in Fig. 2. For the pure DFT calculation 
(f/=0), it shows in Fig. 2(a) that the LDA overbinds the 
compound and underestimates with respect to the ex- 
periment the lattice parameter by ~2%, while the GGA 
calculation give a slight overestimate of ao. After turn- 
ing on the Hubbard U, one can see from Fig. 2(a) that 
for the LDA+t/ approach, although the lattice param- 
eter is still underestimated in a wide range of U , the 
calculated ao for Pu02 improves upon the pure LDA by 
steadily increasing its amplitude with U. In fact, at a 
typical value of [7=4 eV, the LDA+C7 gives ao=5.36 A, 
which is very close to the experiment. On the other hand, 
with increasing U, the underbind effect brought about 
by the GGA+U is somewhat enlarged. As a compari- 
son, at f7=4 eV, the GGA+U gives a =5.47 A, which 
overestimates the experimental data by ~1.3%. Overall, 
both the LDA+L7 and GGA+U results of the lattice pa- 
rameter for the PUO2 AFM phase are comparable with 
experiment by tuning in the calculations the Hubbard U 
around 4 eV. We have also calculated the equilibrium lat- 
tice parameter for the FM and nonmagnetic phases for 
Pu02- The tendency of ao with U for these two phases is 
similar to that for the present AFM phase. For the cal- 
culated bulk modulus Bg of the Pu02 AFM phase, one 
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FIG. 2: Dependence of the lattice parameter (a) and bulk 
modulus (b) of PuCh on U. 



can see from Fig. 2(b) that its value varies with U over 
a rather broad range of 175 to 195 GPa for the GGA+U 
and 200 to 230 GPa for the LDA+U. The LDA result of 
Bq is always higher than the GGA result, which is due 
to the above-mentioned "overbind" characteristics of the 
LDA approach. For the measurements of the equilibrium 
bulk modulus, there are no consistent results to date for 
the AFM Pu02- Here we compared our calculation to the 
experimental result of £?o = 178 GPa reported in Ref. 31 . 
One can see from Fig. 2(b) that the discrepancy between 
the present calculation and the experiment is most dis- 
tinct at U=Q. Both the LDA and the GGA give an over- 
estimate, with the latter more close to the experimental 
data. By turning on the effective Coulomb interaction, 
the amplitude of Bq begins to decrease. At a typical 
value of [7=4 eV, the LDA+C7 gives B =208 GPa while 
the GGA+U gives B =184 GPa. We notice that the re- 
cent hybrid density-funtional calculations 14 predict the 
bulk modulus to be 220 GPa for the antiferromagnetic 
Pu02, comparable with the present pure LDA results. 
For lattice parameter ao it was predicted to be 5.46 in 
RefJ^. To conclude (Fig. 2), comparing with the ex- 
perimental data and the recent hybrid density-functional 
results, the accuracy of our atomic-structure prediction 
for the antiferromagnetic Pu02 is quite satisfactory by 
tuning the effective Hubbard parameter U in a range 3-4 
eV within the LDA/GGA+U approaches. 

Besides the prominent changes in the atomic-structure 
parameters, the most dramatic improvement brought by 
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FIG. 3: The total DOS for the Pu02 antiferromagnetic phase 
computed in the (a) LDA+C7 and (b) GGA+U formalism with 
four selective values of U. The projected DOS for the Pu 5/ 
and O 2p orbitals are also shown. The Fermi level was set to 
be zero. 



the LDA/GGA+[7 when compared to the pure ones is 
in the description of electronic-structure properties. For 
this we have investigated the band structures of the Pu02 
AFM phase with the aim at seeing the fundamental influ- 
ence by the inclusion of the on-site Coulomb interaction. 
The resultant total density of states (DOS) for four se- 
lective values of U are plotted in left (LDA+U) and right 
(GGA+J7) panels in Fig. 3. For more clear illustration, 
the projected DOS for the Pu 5/ and O 2p orbitals are 
also shown in Fig. 3. The Fermi energy Ep has been set 
to be zero. Without accounting for the on-site Coulomb 
repulsion (£7=0), one can see that both two pure DFT 
methods predict an incorrect metallic ground state by 
non-zero occupation of Pu 5/ states at Ep. When switch- 
ing on U, as shown in Fig. 3, the Pu 5/ band begins to 
split at Ep and tends to open a gap A. The amplitude of 
this insulating gap increases with increasing U, see Fig. 
4. Overall the LDA+C7 and GGA+J7 give an equivalent 
description of the one-electron behaviors in a wide range 
of U. At a typical value of [7=4 eV, one can see from Fig. 
3 that the occupied DOS is featured by two well-resolved 
peaks. The narrow one near —2.0 eV is principally Pu 5/ 
in character, while the broad one near —4.0 eV is mostly 
O 2p. These two pronounced peaks have been observed 
in the recent photocmission measurement s] 10 ' 12 . In ad- 
dition, by increasing the amplitude of U , one prominent 
feature occurred in Fig. 3 is the increasing hybridiza- 
tion between Pu 5/ and O 2p occupied states. This in- 
teresting mixing effect disappears in the cases of PU2O3 
(see Fig. 6 below) and UO^, for which the Pu (U) 5/ 
and O 2p occupied bands are well separated. The pres- 
ence of Pu(5/)-Pu(2p) hybridization in Pu02 implies a 
more covalent and stronger metal-ligand mixing than in 
PU2O3 and U02- This phenomenon appears surprising, 
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FIG. 4: The insulating band gap of the PuCh antiferromag- 
netic phase as a function of U for the LDA (filled squares) 
and the GGA (hollow circles). 



given the smaller overlap anticipated in Pu because of 
the smaller radius of the Pu 5/ orbital. Experimentally, 
Butterfield et al. 10 and Gouder et al— have reported 
the thin-film photoemission data for P11O2. The present 
overall picture which emerges from the LDA/GGA+/7 
with properly selective Coulomb repulsion appear to be 
in satisfactory agreement with experiment. We have also 
compared our results given in Fig. 3 with the most recent 
calculations by Prodan et al^ based on newly developed 
screened Coulomb hybrid density functional. The agree- 
ment between our LDA/GGA+J7 (with U ~4 eV) results 
and those in Ref. 14 is also apparent. Interestingly, the 
above-mentioned orbital (Pu 5/ and O 2p) mixing effect 
in Pu02 has also been theoretically predicted by Prodan 
et a/i 14 ' 15 , who hypotheses that the expected stabiliza- 
tion of the Pu 5/ orbital energy relative to U 5/ leads 
to an "accidental" degeneracy between the Pu 5/ and 
O 2p levels, which in the first-order perturbation the- 
ory results in a higher degree of covalency regardless of 
small radius of the Pu 5/ orbital. Therefore, although the 
pure LDA and GGA fail to depict the electronic struc- 
ture, especially the insulating nature and the occupied- 
state character of Pu02, our present results show that by 
tuning the effective Hubbard parameter in a reasonable 
range, the LDA/GGA+J7 approaches will prominently 
improve upon the pure LDA/GGA calculations and thus 
can provide a satisfactory qualitative electronic structure 
description comparable with experiments and the hybrid 
DFT calculation. By further increasing U to 6 eV, one 
can see that the peak near —2.0 eV becomes weak and 
is mostly O 2p, while the peak near —4.0 eV becomes 
stronger and consists equally of Pu 5/ and O 2p orbital. 
This picture of DOS is no longer valid since the peak near 
—2.0 eV has been confirmed to be due to the Pu 5/ con- 
tribution. Thus the LDA/GGA+U approaches with U 
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FIG. 5: Dependence of the equilibrium lattice parameter ao 
(a) and the volume Vb of unit cell (b) of PU2O3 on U. 



as large as 6 eV fails to describe the electronic structure 
of Pu0 2 . 



B. Atomic and electronic structure of PU2O3 

PU2O3 is an insulating oxide of the hexagonal /?- 
type (P3ml) [Fig. 1(b)] with space group no. 164, 
the only phase of the sesquioxide that has been pre- 
pared with stoichiometric composition. Both magnetic 
susceptibility^ and neutron diffraction^ 4 - measurements 
have found PU2O3 to have an AFM structure at tem- 
peratures below 4.2 K, with the Pu moments /x confined 

along the z axis in a simple H 1 — alternation of spins. 

As with PUO2, we have considered the FM, AFM, and 
nonmagnetic phases and then determined the ground- 
state phase by comparing the equilibrium total energies 
of these three phases. At [7=0, the calculated ground 
state is as for Pu02 an incorrect FM metal. By increas- 
ing the amplitude of U, our LDA/GGA+C/ approaches 
correctly predicted the /3-PU2O3 to be in an AFM insu- 
lating phase. The FM-AFM energy crossing occurs at a 
small U of ~1.5 eV. We report in what follows on the 
Pu 2 3 AFM phase. 

The calculated equilibrium lattice parameter ao of 
PU2O3 is plotted in Fig. 5(a) as a function of U. It 
reveals that the relation between ao and U does not fol- 
low a simple monotonic function. The turning point is at 
U =3 eV, below which ao goes up rapidly with U. After 
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FIG. 6: The total DOS for the PU2O3 antiferromagnetic phase 
computed in the (a) LDA+U and (b) GGA+(7 formalism with 
four selective values of U. The projected DOS for the Pu 5/ 
and O 2p orbitals are also shown. 



crossing this turning point, the increase of ao begins to 
slow down with U . Thus the curvature of ao for small 
values of U is more significant than for large values of 
U. The decrease in curvature at large U corresponds to 
the separation of the occupied Pu 5/ band from the un- 
occupied part, i.e., the transition from a metallic to an 
insulating ground state of Pu 2 3 (see below). This fea- 
ture in the increase of ao as a function of U is almost 
the same for the LDA and the GGA. The experimental 
data^ 3 .^ of a =3.841 Ais well fitted at U=i eV for the 
GGA, while the LDA always slightly underestimates ao- 
Another feature shown in Fig. 5(a) is that at small values 
of U below 4 eV the GGA underestimates ao, which is 
contrary to the general experience that in most cases (as 
shown for PUO2), the GGA often gives a slight overesti- 
mate of lattice parameter. This rarely-occurred feature 
may be due to the appearance of the other lattice param- 
eter in /3-PU2O3, i.e., the rario Co/ao for the hexagonal 
crystalline structure. The equilibrium volume Vq of of 
the PU2O3 unit cell (including 5 atoms) as a function of 
U is plotted in Fig. 5(b). The experimental resul t 33 i 34 
of Vq is also given for comparison. Although the ten- 
dency of Vo with U is remarkably opposite for the two 
DFT+f/ methods, the results mostly overlap at a typical 
value of [7=4 eV, at which insulating gap for the PU2O3 
is well formed. The different tendency of Vq with respect 
to U for the LDA and GGA may come from sensitivity of 
the anisotropy in Pu 5/ orbitals to the treatment of the 
exchange-correlation potential. Combining Fig. 5(a) and 
(b) it is expected that both the LDA and the GGA may 
give a satisfactory prediction of the ground-state atomic 
structure for the PU2O3 by tuning U to be near 4 eV. 

The LDA/GGA+J7 total DOS for the Pu 2 3 AFM 
phase are shown in Fig. 6 for four selective values of 
U. The projected DOS for the Pu 5/ and O 2p orbitals 



are also plotted. Both the LDA and GGA predict an in- 
correct metallic ground state for PU2O3 at U=0 by the 
presence of non-zero occupation of Pu 5/ state at the 
Fermi energy Ep. When turning on the on-site Coulomb 
repulsion, the Pu 5/ band begins to split and form an 
insulating gap A at a critical value U=l eV. The gap 
A becomes large with increasing U, as shown in Fig. 
7, from which one can see that the amplitude of A for 
PU2O3 is almost equivalent to that for Pu02 at low U. 
At a typical value of U =4 eV, it reveals in Fig. 6 that 
the occupied DOS is featured by two peaks. The nar- 
row one near —1.5 eV is principally Pu 5/ in character, 
while the broad one around —4.0 eV is mostly O 2p. 
It is encouraging that these two pronounced peaks, as 
well as the overall appearance of the total DOS spec- 
trum, fit well in recent photoemission experiment al 12 
on PU2O3. We have also compared our results given in 
Fig. 6 with the recent calculations by Prodan et alJ^- 
using the hybrid density functional. Our LDA/GGA+C7 
results (with U ~4 eV) for the Pu 2 3 AFM phase are 
in excellent agreement with those in ReL 1 ^. Unlike in 
PUO2, the Pu 5/ and O 2p states in PU2O3 are well sep- 
arated in the DOS spectrum. This feature is similar to 
that of UOz^, which also exhibits two distinct peaks of 
U 5/ and O 2p parentage. Remarkably, the similar trend 
has also been theoretically reported on PU2O3 in Ref. 14 
within the hybrid-dcnsity-functional framework. A con- 
sistent explanation with the Pu(5/)-0(2p) hybridization 
in PUO2 may sustain by understanding the orbital sep- 
aration in PU2O3 as a consequence of the more weakly 
bound Pu 5/ site energy associated with the less highly 
charged Pu 3+ ioni 3 -. With further increasing the effective 
intratomic Coulomb interaction to U =6 eV, as shown in 
Fig. 6, the separation of the Pu 5/ from O 2p projected 
DOS is blurred by the increasing spectrum weight of the 
former around —4 eV, which overlaps largely with the 
O 2p. This no longer accord with the experiment o 10 ' 12 . 
Therefore, as with Pu02, the LDA/GGA+U approaches 
with U as large as 6 eV fail to describe the electronic 
structure of PU2O3. 



C. Oxidation reaction energy 



Oxidation of PU2O3 via the reaction 



1 



Pu 2 3 + -0 2 -► 2Pu0 2 



(2) 



leads to formation of stoichiometric Pu02- The depen- 
dence of the transformation reaction energy on U is pre- 
sented in Fig. 8. One can see that both the LDA and the 
GGA show the same dependence of reaction energy on 
the on-site Coulomb interaction. That is, at small val- 
ues of U which correspond to the metallic ground state 
for both Pu02 and PU2O3, the reaction energy is inde- 
pendent of U. Above the metallic-insulating transition, 
our calculated reaction energy decreases linearly with in- 
creasing U. The reason for this behavior is that a high 
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FIG. 7: The insulating band gap of the PU2O3 antiferromag- 
netic phase as a function of U for the LDA (filled squares) 
and the GGA (hollow circles). 
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U favors localization and thus facilitates the transition. 
Density functional theory is known to overestimate the 
binding energy of O2 and this should result in an un- 
derestimation of the present reaction energy via the Eo 2 
term. Consequently we cannot expect a perfect agree- 
ment with experiments for the present reaction energy. 



However, this error is independent of any conditions in 
the plutonium oxide and thus can be remedied by shifting 
the energy of O2 so as to give the experimental binding 
energy. In the LDA the O2 binding energy is overesti- 
mated by 1.2 eV/0.5 O2 and in the GGA the correspond- 
ing number is 0.8 eV/0.5 O2. The GGA always predicts 
a lower value of the reaction energy, as seen from Fig. 9. 
IV. CONCLUSIONS 

We have studied the structural, electronic, and thermo- 
dynamic properties of the antiferromagnetic PUO2 and 
PU2O3 within the LDA+J7 and GGA+C7 frameworks. 
The atomic structure, including lattice parameters and 
bulk modulus, and the one-electron behaviors of these 
kinds of plutonium oxides have been systematically in- 
vestigated as a function of the effective on-site Coulomb 
repulsion parameter U. We find that both the LDA+C7 
and GGA+C7 considerably improves upon the traditional 
density functionals, providing a first-principles descrip- 
tion of plutonium oxides in satisfactory qualitative agree- 
ment with experiment. Also our present results are well 
comparable to those obtained through newly developed 
hybrid DFT method. Specially, from the LDA/GGA+C/ 
study of the lattice parameter of Pu02 we find that the 
experimental data of ao can be gradually approached 
by steadly increasing U to be in an acceptable range 
around 4 eV. The incorrect metallic ground state at 
purely LDA or GGA (U=0) for both Pu0 2 and Pu 2 3 
can be readily corrected by a systematic inclusion of non- 
zero U, which forces the Pu 5/ band to split at the 
Fermi level and thus drives the metallic-insulating tran- 
sition. The insulating band gaps for PuC>2 and PU2O3 
have been shown as a function of U. The oxidation re- 
action PU2O3+O.5O2 — >2PuC>2 has also been studied by 
systematically calculating the reaction energy as a func- 
tion of U. Our results show that the oxidation process 
of the PU2O3 is an exothermic reaction, which is mostly 
responsible for the experimentally observed^ plutonium 
pyrophoricity at 150°C— 200°C. Also we have shown that 
above the metallic-insulating transition, the reaction en- 
ergy decreases with increasing U for the LDA and the 
GGA schemes. We expect these calculated results are 
useful for the future studies on the surface oxidation and 
corrosion of metallic plutonium. 
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